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Summary. A new generation of observatories is looking for gravitational waves. 
These waves, emitted by highly relativistic systems, will open a new window for ob- 
servation of the cosmos when they are detected. Among the most promising sources 
of gravitational waves for these observatories are compact binaries in the final min- 
utes before coalescence. In this article, we review in brief interferometric searches for 
gravitational waves emitted by neutron star binaries, including the theory, instru- 
mentation and methods. No detections have been made to date. However, the best 
direct observational limits on coalescence rates have been set, and instrumentation 
and analysis methods continue to be refined toward the ultimate goal of defining 
the new field of gravitational wave astronomy. 



1 Introduction 

It is no great exaggeration to say that the advent of a new type of astronomy 
is imminent. Gravitational wave astronomy is predicated on the observation 
of the cosmos not with a new band of the electromagnetic spectrum, but 
rather via a whole new spectrum, the spectrum of gravitational waves. As 
such it has the potential to revolutionize our understanding of the Universe 
because it will allow us to access phenomena which are electromagnetically 
dark or obscured [1]. The impetus for this revolution is the recent construction 
and operation of a new generation of gravitational wave detectors based on 
interferometry [2]. Within a decade, these interferometers will be reaching 
sensitivities at which gravitational wave observations should become routine. 

There are four categories of gravitational wave signals which ground-based 
interferometers are currently trying to detect: quasi-periodic signals, such as 
those expected from pulsars [3, 4, 5, 6, 7], stochastic background signals, such 
as remnant gravitational waves from the Big Bang [8, 9, 10, 11, 12], unmodeled 
burst signals, such as those that might be emitted by supernovae [13, 14, 
15, 16, 17, 18, 19], and inspiral signals, such as those from neutron star or 
black hole binaries [20, 21, 22, 23, 24, 25, 26]. In this article, we will only 
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be concerned with the last of these searches, and in particular the search for 
neutron star binaries [20, 23, 25, 26]. 

This article is organized as follows: Since many astronomers may not be 
very familiar with gravitational waves and the effort to use them for astron- 
omy, Sect. 2 is devoted to background material. This includes a simple de- 
scription of what gravitational waves are and how they are generated, a brief 
history of the instruments and efforts to detect them, and some background 
on the relevant aspects of neutron star binaries and the gravitational waves 
they produce. Section 3 describes in some detail how searches for gravitational 
waves from neutron star binaries have been conducted, including a description 
of the data, the data analysis methods employed, and coincidence vetoes. 

Section 4 reviews the published upper limits that have been placed on the 
rate of neutron star binary coalescence in our galactic neighborhood by inter- 
ferometric detectors. Also included is a description of the statistical analysis 
used to place these upper limits. In Sect. 5, we discuss prospects for bet- 
ter upper limits and discuss some of the astrophysics that might be done 
by interferometric gravitational wave detectors when they reach a sensitivity 
where gravitational waves from neutron star binaries are regularly observed. 
Concluding remarks are found in Sect. 6. 

2 Background 

2.1 Gravitational Waves 

One of the many remarkable predictions of Einstein's general theory of rel- 
ativity is the existence of gravitational waves (GWs) [27]. Einstein himself 
elucidated the theoretical existence of GWs as early as 1918 [28]. Today, how- 
ever, GWs have still not been directly measured, although the measurements 
of the binary pulsar PSR1913+16 [29, 30, 31], (discovered by Hulse and Taylor 
and for which they won the Nobel prize) leave little doubt that GWs do, in 
fact, exist. 

The fundamental factor that has led to our failure to directly measure GWs 
is the exceptional weakness of the gravitational coupling constant. This causes 
GWs to be too feeble to detect unless produced under extreme conditions. In 
particular, gravitational waves are produced by accelerating masses (a more 
exact formulation of this statement appears in Sect. 2.3). Thus, for the highest 
GW amplitudes, we seek sources with the highest possible accelerations and 
masses. As a result, astrophysical objects are the most plausible sources of 
GWs 3 . 

3 The exception to this statement is the big-bang itself, which should lead to a 
stochastic cosmological background of GWs, as mentioned in Sect. 1. Gravita- 
tional wave observations have just begun to bound previously viable theoretical 
models of this background [10]. We will not be considering this background further 
in this review. 
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Further restrictions on viable sources are imposed by our detectors. For 
instance, for Earth based instruments, even with the best seismic isolation 
technologies currently available [32, 33], noise from seismic vibrations lim- 
its large-scale precision measurements to frequencies above ^30 Hz. Through 
causality, this time-scale limitation implies a maximum length-scale at the 
source of 10 4 km. Given the Chandrasekhar limit on the mass of a white 
dwarf star [34, 35] and the white dwarf mass-radius relationship this is ap- 
proximately the minimum length-scale for white dwarf stars [36]. 

If searches are restricted to objects more compact than white dwarfs, then 
within the bounds of current knowledge, gravitational wave astronomy with 
ground-based interferometers is limited to black holes and neutron stars as 
sources. According to our current understanding of astrophysical populations, 
these objects are relatively rare. Thus, the probability of finding them in our 
immediate stellar neighborhood are small, and any realistic search must be 
sensitive to these sources out to extragalactic distances to have a reasonable 
chance of seeing them in an observation time measured in years. Is this a 
reasonable prospect? The answer is yes, but to understand why, it behooves 
us to first understand a bit more about what GWs are and how they might 
be measured. 

Gravitational waves arriving at Earth are perturbations of the geometry 
of space-time. Heuristically, they can be understood as fluctuations in the 
distances between points in space. Mathematically, they are modeled as a 
metric tensor perturbation h a p on the flat spacetime background. Linearizing 
the Einstein field equations of general relativity in h a p, we find that this 
symmetric four-by-four matrix satisfies the the wave equation, 



in an appropriate gravitational gauge. Here V 2 is the usual Laplacian oper- 
ator, G is the gravitational constant, and T a p is the stress energy tensor, 
another symmetric four-by-four matrix which encodes information about the 
energy and matter content of the spacetime. 

From (1), it is obvious that GWs travel at c, the speed of light. To under- 
stand the production of gravitational waves by a source, we solve (1) with the 
stress-energy (T a p) of that source on the right-hand-side. We will discuss this 
in more detail in Sect. 2.3. First, however, we wish to consider the propagation 
of gravitational waves. 

For the propagation of gravitational waves, we require solutions to the ho- 
mogeneous (T a p = 0) version of (1). As usual, such solutions can be expressed 
as linear combinations of the complex exponential functions 




(1) 



h' 



p = A a pexp(±ik fl x f *). 



(2) 



Here, A a p is a matrix of constant amplitudes, fc M = (— u>, k) is a four vector 
which plays the role of a wave vector in four dimensions, and = (t, x) arc 
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the spacetime coordinates. Above, and in what follows, we use the Einstein 
summation convention that repeated indices, such as the /i in k^x^ = —cut + 
k x, indicate an implicit summation. It can be shown that for gauges in which 
(1) holds that 

k^p = 0. (3) 

In words, this means that the wave vector is orthogonal to the directions in 
which the GW distorts spacetime, i.e. the wave is transverse. 

Since h a p is a four-by-four matrix, it has 16 components. However, because 
it is symmetric, only 10 of those components are independent. Further, (3) 
imposes four constraints on h a p, reducing the number of free components to 
six. One can use remaining gauge freedom to impose four more conditions 
on h a p. There are therefore only two independent components of the matrix 
h a p. Details can be found in any elementary textbook on General Relativity, 
such as [37]. 

The two independent components of h a p are traditionally called h+ and 
h x . These names are taken from the effect that the components have on a ring 
of freely moving particles laying in the plane perpendicular to the direction 
of wave propagation, as illustrated in Fig. 1. The change in distance between 




Fig. 1. Distortion of a ring of freely falling dust as a gravitational wave passes 
through. The wave is propagating into the page. From left to right are a series of 
four snapshots of the distortion of the ring. The top row are distortions due to 
h+. The bottom row are distortions due to hx- The snapshots are taken at times 
t = 0, t = T/4, t = T/2 and t = 3T/4 respectively, where T is the period of 
the gravitational wave. The relative phase between h + and h x corresponds to a 
circularly polarized gravitational wave. 



particles, 6£, is proportional to the original distance between them, £, and 
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the amplitude of the gravitational waves, A a p. For a point source, which all 
astrophysical sources will effectively be, the amplitude decreases linearly with 
the distance from the source. For astrophysical source populations from which 
gravitational wave emission have been estimated, the typical gravitational 
wave strain, h ~ 2S£/£, at a detector at Earth would be expected to be less 
than or of the order of 10~ 21 [1]. Through interferometry, it is possible to 
measure 81 <~ 1CP 18 m. Thus, interferometers of kilometer scales are required 
to have any chance of measuring these sources. 

It might seem that the challenge of attempting to measure gravitational 
waves is so daunting as to call into question whether it is worthwhile at all. 
However, there are several factors which make the measurement of gravita- 
tional waves attractive. First, astrophysical gravitational wave sources include 
systems, such as black hole binaries, which are electromagnetically dark. Grav- 
itational waves may therefore be the best way to study such sources. Second, 
since gravitation couples weakly to matter, gravitational waves propagate es- 
sentially without loss or distortion from their source to the detector. Thus, 
sources obscured by dust or other electromagnetically opaque media may still 
be observed with gravitational waves. Also, interferometers behave as ampli- 
tude sensing devices for GWs (like antennas), rather than energy gathering 
devices (like telescopes), leading to a 1/r fall-off with distance, rather than 
the more usual 1/r 2 fall-off [2]. 

But perhaps the most compelling reason for pursuing the measurement 
of gravitational waves is that they constitute an entirely new medium for 
astronomical investigation. History has demonstrated that every time a new 
band of the electromagnetic spectrum has become available to astronomers, 
it has revolutionized our understanding of the cosmos. What wonders, then, 
await us when we start to see the Universe through the lens of gravitational 
waves, (which will surely begin to happen within the next decade as GW 
detectors continue their inevitable march toward higher sensitivities)? Only 
time will tell, but there is every reason to be optimistic. 

In this article, we concentrate on one of the many sources of gravitational 
waves for which searches are ongoing - binary neutron star (BNS) systems. 
In particular, the sensitive frequency band of ground-based interferometers, 
which is approximately 40 Hz to 400 Hz, dictates that we should be inter- 
ested in neutron star binaries within a few minutes of coalescence [38] . These 
sources hold a privileged place in the menagerie of gravitational waves sources 
that interferometers are searching for. The Post-Newtonian expansion, a gen- 
eral relativistic approximation method which describes their motion, gives 
us expected waveforms to high accuracy [39, 40, 41]. They are one of the few 
sources for which such accurate waveforms currently exist, and they are there- 
fore amenable to the most sensitive search algorithms available. Furthermore, 
while the population of neutron star binaries is not well understood, there are 
at least observations of this source with which to put some constraints on the 
population [42] . These two factor give BNS systems one of the best (if not the 
best) chance of discovery in the near future. 
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2.2 Gravitational Wave Detectors 

Gravitational wave detectors have been in operation for over forty years now. 
However, it is only in the last five years or so that detectors with a non- 
negligible chance of detecting gravitational waves from astrophysical sources 
have been in operation. The first gravitational wave detector to operate was 
built by Joseph Weber [43] . ft consisted of a large cylinder of aluminum, two 
meters long and a meter in diameter, with piezoelectric crystals affixed to 
either end. 

The fundamental idea for such detectors is that if a strong enough gravita- 
tional wave was to pass by, that it would momentarily reduce the interatomic 
distances, essentially compressing the bar, and setting it ringing, like a tuning 
fork. The ringing would create electric voltages in the piezoelectric crystals 
which could then be read off. Of course, as with a tuning fork, the response 
of the apparatus is greatly increased if it is driven at its resonant frequency 
(about 1660Hz, for Weber's bars). 

Weber's bar was isolated from seismic and electromagnetic disturbances 
and housed in a vacuum. He attributed the remaining noise in his instrument 
to thermal motion of the aluminum atoms. This noise limited Weber's mea- 
surements to strains of h^l0~ 16 , about five orders of magnitude less sensitive 
than the level now believed necessary to make the probability of detection 
non-negligible. Nonetheless, by 1969, Weber had constructed two bars and 
had observed coincident events in them although they were housed approx- 
imately fOOO miles apart. He calculated that his noise would create some of 
these events at rates as low as one per thousand years, and subsequently 
published his findings as "good evidence" for gravitational waves [44] . 

Many groups followed up on this and subsequent claims of gravitational 
wave detections made by Weber. No other group was ever able to reproduce 
these observations, and it is now generally agreed that Weber's events were 
spurious. Nonetheless, this launched the era of gravitational wave detectors, as 
resonant mass detectors (as Weber-like bars are now called) began operating 
in countries around the globe. Today, there is a network of these detectors op- 
erated under the general coordination of the International Gravitational Event 
Collaboration [45] . Technical advances have led to considerable improvements 
in sensitivity over the past decades. They continue to be rather narrow band 
detectors, however, and they are not the most sensitive instruments in oper- 
ation today. 

That distinction belongs to intcrfcromctric gravitational wave detectors [2] . 
The basic components of these detectors are a laser, a beam splitter to divide 
the laser light into two coherent beams, hanging mirrors to reflect the laser 
beams, and a light sensing diode, as shown in Fig. 2. The mirrors act essen- 
tially as freely moving particles in the horizontal directions. If a sufficiently 
strong gravitational wave with a vertical wave-vector component impinges, it 
shortens the distance between the beam splitter and one of the mirrors, and 
lengthens the distance between the beam splitter and the other mirror. This 
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Fig. 2. Schematic diagram of an interferometric gravitational wave detector. The 
beam splitter is coated to allow half the light to be transmitted to one of the mirrors, 
and the other to be reflected to the other mirror. Real interferometric gravitational 
wave detectors are much more sophisticated, including frequency stabilization of 
the laser, a second mirror on each arm between the beam splitter and the end 
mirror to create Fabry-Perot cavities, and control feedback loops which "lock" the 
interferometer onto an interference fringe. Thus, rather than measuring the current 
from the photo-diode directly, the gravitational wave signals are encoded in the 
feedback loop voltages needed to maintain the lock of the interferometer. These and 
many other enhancements are necessary to reach the required sensitivity level. 

is registered as a shifting interference pattern by the light sensing diode, thus 
detecting the gravitational wave. 

More specifically, what is measured is a quantity proportional to the strain 
on the detector, s :— (Sx — Sy)/£, where Sx and Sy are the length changes 
in the two equal length arms of the interferometer, traditionally called the x 
and y arms, and £ is the unperturbed length of each arm. If there were no 
noise, then in the special case that h + or h x was aligned with the arms, we 
would have Sx = —Sy = SI and the measured quantity would be proportional 
to 2S£/£. For a more general alignment, but still in the absence of noise, some 
linear combination of h + and h x is measured 



where F+ and F x are called beam pattern functions. They project the gravi- 
tational wave components, h + and h Xl onto a coordinate system defined by 
the detector, and are functions of the Euler angles (0, <j>, tp) which relate this 
coordinate system to coordinates which are aligned with the propagating GW. 

In the case where there is noise, Sx and Sy are sums of the gravitational 
wave displacements and the noise displacements, so that the detector strain 
is 



h(t) := F + (e,<t>,4>)h+(t)+F x (9,<t>,i))h x (t), 



(4) 
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s(t) = h(t) + n(t) (5) 

where n(t) is the noise contribution. If the noise component can be kept 
from dominating the signal component, there is a reasonable chance that 
the gravitational wave can be detected. For a more complete and detailed 
account of interferometric detection of gravitational waves, as well as many 
other aspects of gravitational wave physics, we recommend [1]. 

The idea of using an interferometer as a gravitational wave detectors was 
explored repeatedly but independently by several different researchers over 
approximately 15 years [46, 47, 48, 49]. The first prototype of an interfer- 
ometric gravitational wave detector was built by one of Weber's students, 
Robert Forward, and collaborators in 1971 [50, 51]. The advantages of this 
idea were immediately understood, but, as mentioned above, an understand- 
ing also emerged that kilometer-scale interferometry would be needed. Thus, 
these detectors needed to be funded, built and operated as coordinated efforts 
at the national or international level. 

To date, there have been six large-scale (100 m plus) interferometers oper- 
ated at five sites. Three of the these are located in the United States and con- 
stitute the Laser Interferometer Gravitational-wave Observatory (LIGO) [52, 
53] . There is one four kilometer instrument at each of the LIGO-Hanford site 
in Washington State [54] and LIGO-Livingston site in Louisiana [55] (HI and 
LI respectively), and a two kilometer instrument, housed in the same en- 
closure, at the LIGO-Hanford Observatory (H2). The other interferometers 
arc the three kilometer Virgo instrument, built in Italy by a French-Italian 
collaboration [56], the 600 meter German-English Observatory (GEO600) in 
Germany [57], and the TAMA observatory, a 300 meter instrument located in 
and funded by Japan [58]. 

Other large-scale interferometers are in various stages of planning, al- 
though funding has not been secured for them. One of the most interesting 
is LISA, a planned joint NASA-ESA mission, which would consist of two in- 
dependent million-kilometer-scale interferometers created by three satellites 
in solar orbit [59, 60]. The larger scale and freedom from seismic noise will 
make this instrument sensitive at much lower frequencies than its Earth-based 
brethren. 

Because they are currently the most sensitive gravitational wave detectors 
in the world, and the ones with which we are most familiar, this article will 
often use the particular example of LIGO and its methods to illustrate our 
discussion. LIGO began construction in 1994 and was commissioned in 1999. It 
began taking scientific data on 23 August, 2002. That data taking run, called 
Science Run One (or SI for short), ended 9 Sept., 2002. There have been four 
subsequent science runs: S2 from 14 Feb. to 14 Apr. 2003, S3 from 31 Oct. 
2003 to 9 Jan. 2004, S4 from 22 Feb. to 23 Mar. 2005, and finally S5 has been 
ongoing from Nov. 2005 and is expected to end in fall, 2007. At the conclusion 
of S5, the LIGO interferometers are scheduled for component enhancements 
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which are expected to double the sensitivity of the instrument [61, 62]. The 
upgrade process is scheduled to last approximately a year. 

As with all the new interferometers, initial runs were short and periods 
between them were long because scientists and engineers were still identifying 
and eliminating technical and environmental noise sources which kept the in- 
struments from running at their design sensitivity. As designed, interferometer 
sensitivity is expected to be bounded by three fundamental noise sources [49] . 
Below approximately 40 Hz, seismic noise transmitted to the mirrors through 
housing and suspension dominates. Between approximately 40 and 200 Hz, 
thermal vibrations in the suspension system for the mirrors dominates. Fi- 
nally, above approximately 200 Hz, the dominant contribution comes from 
photon shot noise associated with counting statistics of the photons at the 
photodiode. These three fundamental noise regimes, which contribute to the 
"noise floor" of a detector, are described in the caption of Fig. 3. 



Best Strain Sensitivities for the LIGO Interferometers 
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Fig. 3. Interferometer noise in each of the five LIGO science runs. The solid thick 
black curve is design goal of LIGO. There are two changes of slope indicated by 
dashed thick black lines. These correspond to changes in the dominant noise source 
- the three noise regimes are marked. Over four years of commissioning, the noise 
floor was reduced by approximately three orders of magnitude. 
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Apparent in Fig. 3 is the remarkable progress that has been made in low- 
ering the noise floor and the other noises superimposed upon it in going from 
SI to S5. Indeed, at S5, LIGO matches or exceeds its design specification over 
most frequency bands, including the most sensitive region from approximately 
100-300 Hz. When S5 is over and S5 data is analyzed, it will provide by far 
the best opportunity to date to detect gravitational waves. 

The degree to which improvements in the noise floor such as those in Fig. 3 
correspond to increase detection probabilities depend on the shape of the noise 
floor and the signal being sought. For neutron star binary inspiral signals, one 
can devise a figure of merit called the inspiral range. This is the distance, 
averaged over all sky positions and orientations, to which an instrument can 
expect to see a 1.4-1.4 M Q binary with signal-to-noise ratio (defined below) of 
8. For S2, typical inspiral ranges for LI were just below 1 Mpc. For the LIGO 
design curve, the inspiral range is approximately 14 Mpc [63]. 

2.3 Neutron Star Binaries as Gravitational Wave Emitters 

In the preceding subsections, we discussed the propagation and detection of 
gravitational waves. It is perhaps useful now to say a few words about the 
production of gravitational waves before diving into the particular sources of 
interest for this article, binary neutron star systems. 

We have discussed solutions to the homogeneous (T a p = 0) version of (1), 
which describe the propagation of GWs in the far wave zone. The generation of 
gravitational waves requires a source, and is therefore described by (1) with 
T a p^0. In fact, to lowest order, the relevant property of the source is the 
mass density of the source, p(t,x), or more specifically, the mass quadrupole 
moments. 

The relevant quadrupole moments depend on the direction from the source 
at which the gravitational wave is detected. In the limit of a negligibly gravi- 
tating source, the moments are given by the integrals 



Here, x — {x, y, z} are Cartesian coordinates centered on the source with the 
2-axis being defined by the direction to the detector and t is time. In terms of 
these integrals, the relevant components of the solution to (1) with source is 




(6) 



(7) 



M*) = - ? ^ + (*-r). 



(8) 



and likewise for h x (t). Here r is the distance between the source and the 
detector. 
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While higher order multipole moments of the mass distribution can con- 
tribute to the radiation, for most systems the quadrupole will dominate. Fur- 
ther, the mass monopole and dipole moment will not contribute any gravita- 
tional waves. Thus, such events as a spherically symmetric gravitational col- 
lapse and axially symmetric rotation do not emit any gravitational radiation. 
On the other hand, a rotating dumbell is an excellent emitter of gravitational 
waves, making binary systems potentially amongst the brightest emitters of 
gravitational waves in the Universe. 

In general, the better the information about the gravitational wave signal, 
the better will be our chances of detecting it. This, in turn, prods us to find 
the best possible model of the dynamics of a GW emitting system. For the 
gravitational wave sources suitable for Earth-based interferometric detectors, 
this proves to be difficult in general. For instance, the physics of core-collapse 
supernovae is not understood at the level of detail needed to accurately predict 
gravitational waveforms. This is also true of colliding black holes or neutron 
stars, which should emit large amounts of energy in gravitational waves. In 
fact, there are only two sources for ground-based interferometers, inspiralling 
binary neutron stars and perturbed intermediate mass (~100 M ) black holes, 
for which the physics is well enough understood that it is believed that the 
waveforms calculated provide a high degree of confidence for detection 4 . 

Of these two potential sources, we at least know with great certainty that 
neutron star binaries, like the Hulse- Taylor binary, exist [42, 64, 65]. Accord- 
ing to current thinking, these binaries are formed from the individual collapse 
of binaries of main sequence stars [66] . If one can establish a population model 
for our galaxy, therefore, one should be able to deduce the population outside 
the galaxy by comparing the rate of star formation in our galaxy to elsewhere. 
One measure of star formation is blue light luminosity of galaxies (corrected 
for dust extinction and reddening) [67, 68]. Thus, once a galactic population 
model is established, given a noise curve for an instrument it is possible to es- 
timate the rate at which gravitational waves from binary neutrons star (BNS) 
systems out to extragalactic distances will be seen. 

To develop a population model, astrophysicists use Monte Carlo codes 
to model the evolution of stellar binary systems [69, 70]. To cope with un- 
certainties in the evolutionary physics of these systems, a great many (~30) 
parameters must be introduced into the models, some of which can cause the 
predicted rates to vary by as much as two orders of magnitude. This can be 
reduced somewhat by feeding what is known about BNS systems from ob- 
servation into the models [71]. However, only seven BNS systems have been 

4 To be fair, although the source mechanics of rotating neutron stars, like pulsars, 
are not well understood, the gravitational waves they produce if there is some 
asymmetry about the axis of rotation will be quasi-periodic. We therefore do 
not need to model these sources well to search for their gravitational waves. 
Similarly, although a stochastic background of gravitational waves in intrinsically 
unknowable in detail, a search for these waves can be optimized if their statistical 
properties are known. 
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discovered in the galactic disk. Furthermore, the most relevant for detection 
rates with interferometric GW detectors are the four BNS which are tight 
enough to merge within 10 Gyrs, since, as mentioned above, only those within 
minutes of merger can be observed. Thus, there is not much information to 
feed into these models, and estimates can still vary widely. The current best 
estimate is that LIGO should now be able to observe of the order of one 
BNS system approximately every hundred years, although uncertainties ex- 
tend this from a few per thousand years to almost one every ten years [66]. 
An improvement in LIGO's ability to detect BNS systems by a factor of three 
would therefore raise the most optimistic case to almost one BNS signal per 
year. 

With projected rates this low, it is essential to search for BNS systems with 
the highest possible detection efficiency. Since, in general, the more informa- 
tion that can be fed into the detection algorithm, the more efficient it will be, 
it is important to have a high accuracy theoretical prediction for the wave- 
form. In the case of BNS systems, this prediction is provided by the restricted 
post-Newtonian approximation [39, 40, 41]. The post-Newtonian formalism 
uses an expansion in orbital velocity divided by c, the speed of light. Since 
this is a slow motion approximation, and the orbital velocity of the binary 
constantly increases during inspiral, this approximation becomes worse as the 
binary evolves. However, the accuracy is still good for BNS systems when they 
are in the interferometric detection band. Furthermore, we need only deal with 
circular orbits since initially eccentric orbits are circularized through the GW 
emission process [72]. Finally, one can safely ignore spin terms [73] and finite 
size effects [74]. Figure 4 shows the post-Newtonian prediction for a neutron 
star binary in the sensitive frequency band of an interferometric detector. 

The general features of this waveform are easily understood. From (7), 
it is apparent that the mass quadrupole moment will be proportional to the 
square of the orbital radius, a 2 . The second time derivative of this quadrupole 
moment therefore goes as a 2 w 2 , and from (8) we have that h^a 2 ui 2 . For a Ke- 
plerian orbit, which relativistic orbits approximate when the Post-Newtonian 
approximation is valid, aj 2 ^a~ 3 and therefore h^a^ 1 . Thus, when the binary 
emits GWs, which carry away orbital energy, it experiences a progressive tight- 
ening, leading to an increased frequency and amplitude. If the two neutron 
stars are of equal mass, then the quadrupole configuration repeats itself ev- 
ery twice every orbital period, and the gravitational wave frequency is twice 
the orbital frequency. Because the two masses are comparable in every BNS 
system, the dominant frequency component is always at twice the orbital fre- 
quency. 
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(a) Second-order restricted post-Newtonian waveforms from a 1.4-1.4 M Q neutron 
star binary at 1 Mpc 
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(b) Closeup of above waveforms 



Fig. 4. Second-order restricted post-Newtonian h+ and h x waveforms from an 
optimally located and oriented 1.4—1.4 M© neutron star binary system at 1 Mpc. The 
top panel, (a), shows approximately the portion (with frequency sweeping through ~ 
40-400 Hz) of the waveform visible to a ground-based interferometric GW detector. 
The bottom panel, (b), shows a closeup of the last 0.2s of panel (a). From the 
closeup, we see that the waveforms are sinusoids sweeping upward in both frequency 
and amplitude as the binary companions inspiral toward one another. The two 
phases of the gravitational wave, h + and h x , are 90° out of phase. Thus, this is a 
circularly polarized gravitational wave. 



3 Search Method 

3.1 Interferometric Data 

As was mentioned in the caption of Fig. 2, gravitational wave interferometers 
are not the simple Michelson interferometers portrayed in that figure. Rather, 
they have a number of important and sophisticated refinements to that basic 
configuration, all designed to increase instrument sensitivity and reliability 
[2, 75]. For the purpose of understanding the data, the most important en- 
hancement is that, when the IFOs are operational, their arm- lengths are held 
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fixed, or locked. This is done by measuring the light at the antisymmetric port 
(the diode in Fig. 2) and applying feed-back controls to the mirrors through 
magnetic couplings [76, 77]. By measuring the feed-back loop voltage, it is 
possible to monitor how much the arm-lengths would have changed had the 
mirrors been free. This allows the suppression of pendulum modes which are 
cither spurious or which might damage the optics if driven at resonance. The 
feedback control voltage is therefore the gravitational wave channel of the 
interferometer. 

In order to convert this voltage into detector strain, the frequency de- 
pendent transfer function between these two quantities must be applied. The 
overall shape of this transfer function in the frequency domain is determined 
by a model of the instrument, however, the overall scaling of the transfer func- 
tion must be measured. This can be done, for instance, by driving the mirrors 
at specific frequencies and measuring the response of the control loop volt- 
age [78]. At the frequencies at which these calibration signals are injected, one 
can see sharp line features in the noise spectra of the interferometers. These 
are some of the lines found in Fig. 3. Other lines are caused by resonances in 
the wires suspending the mirrors. Lines at multiples of 60 Hz are caused by 
coupling of the electronic subsystems to the power grid mains. Fortunately, 
for the BNS search, these lines are not very problematic because the signal 
sweeps through frequencies and therefore never corresponds to one of these 
noisy frequencies for long. 

As well as the gravitational wave channel, there are usually a myriad 
of other channels monitoring the physical environment of the interferome- 
ter (seismic, magnetic, etc) and the internal status of the instrument (mirror 
alignment, laser power, etc). These can be used to veto epochs of gravitational 
wave data which are unreliable. Furthermore, if one has access to simultaneous 
data from multiple interferometers (like LIGO does with its three interferom- 
eters, all of which have similar alignment 5 and two of which are co-located), 
most false alarms can be eliminated by requiring coincidence of received sig- 
nals in all interferometers. 

Although noise from an interferometer that is functioning well is primarily 
Gaussian and stationary, there are occasional noise excursions, called glitches. 
The causes of some of these glitches are not apparent in the auxiliary channels. 
These constitute the primary source of background noise in the search for 
signals from neutron star binaries. 

3.2 Matched Filtering and Chi-Squared Veto 

When the expected signal is known in advance and the noise is Gaussian and 
stationary, the optimal linear search algorithm is matched filtering [79] . The 
idea behind matched filtering is to take the signal, and data segments of the 

The detectors in Washington are somewhat misaligned with the detector in 
Louisiana due to the curvature of the Earth. 
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same length as the signal, and treat them as members of a vector space. As 
with any two vectors in a vector space, the degree to which the signal and a 
data vector overlap is calculated using an inner product. 

To be more precise, consider detector strain s(t) and a signal h(t) that 
lasts for a duration of T. If the signal arrives at the detector at time to, then 
the detector strain can be written 

s(t) = f h (t-to) + n(t), ta<t<t + T 
\n(t), otherwise 

where n(t) is the detector noise. For this paragraph, we will assume that, apart 
from being stationary and Gaussian, n(t) is white (same average power at all 
frequencies) for simplicity. Then, the matched filter output, ((t), is given by 

C(t) = 2 / h(r) s(t + t) dr. (10) 
J o 

At time t = to, we have 



T i-T 

2r 



C(*o) = 2/ h z { T )d T + 2 h(T)n(t Q + T)dT. (11) 
Jo Jo 

Let us denote the first and second integrals in (11) by I\ and hito) respec- 
tively. Clearly, the integrand of l\ is deterministic and positive everywhere. 
However, the integrand of li is stochastic. The average of li over all noise real- 
izations vanishes. In other words, on average C(^o) = Ii when there is a signal 
starting at time to- On the other hand, when there is no signal, ((t) = hit). 
Denoting the standard deviation of I 2 over all noise realizations by a, we 
define the signal-to-noise ratio (SNR) for the data to be 

g(t) := \((t)\/a. (12) 

Clearly, at time to the expected value of the SNR is g(to) = h/cr. Thus, if the 
signal is strong enough that I\ is several times larger than cr, there is a high 
statistical confidence that it can be detected. 

In practice, it is preferable to implement the matched filter in the frequency 
domain. Thus, rather than a stretch of data s(t), one analyzes its Fourier 
transform 

/oo 
e- 2 " ift s(t)dt, (13) 
— oo 

where / labels frequencies. This has several advantages: first, it allows for 
the non-white noise spectrum of interferometers (cf. Fig. 3) to be more easily 
handled. Second, it allows the use of the stationary phase approximation to 
the restricted post-Newtonian waveform [1, 80], which is much less computa- 
tionally intensive to calculate, and accurate enough for detection [81]. Third, 
it allows one to easily deal with one of the search parameters, the unknown 
phase at which the signal enters the detector's band. 
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In the frequency domain, the matched filter is complex and takes the form 
z(t) = x(t) + iy(t) = 4 j~ ^£]jf df, (14) 

where S n (f) is the one-sided noise strain power spectral density of the detector 
and the * superscript denotes complex conjugation. It can be shown that the 
variance of the matched filter due to noise is 

2 A P h'(f)Hf) 

° =4 Jo^wr df - (15) 

In terms of z(t), the SNR is given by 

e (t) = \z(t)\/a. (16) 

Note that a and z(t) are both linear in their dependence on the signal template 
h. This means that the SNR is independent of an overall scaling of h(t), which 
in turn means that a single template can be used to search for signals from 
the same source at any distance. Also, a difference of initial phase between 
the signal and the template manifests itself as a change in the complex phase 
of z{t). Thus, the SNR, which depends only on the magnitude of the matched 
filter output, is insensitive to phase differences between the signal and the 
template. 

Equations (14-16) tell us how to look for a signal if we know which signal 
to look for. However, in practice, we wish to look for signals from any neutron 
star binary in the last minutes before coalescence. Because, as mentioned 
above, finite-size effects are irrelevant, a single waveform covers all possible 
equations of state for the neutron stars. Likewise, as stated above, the spinless 
waveform will find binaries of neutron stars with any physically allowable spin. 
Further, as just discussed, a single template covers all source distances and 
initial signal phases. However, a single template does not cover all neutron 
star binaries because it does not cover all masses of neutron stars. 

Population synthesis models for neutron star binaries indicate that masses 
may span a range as large as ~1-3M . Since mass is a continuous parameter, 
it is not possible to search at every possible mass for each of the neutron stars 
in the binaries. However, if a signal is "close enough" to a template, the loss 
of SNR will be small. Thus, by using an appropriate set of templates, called 
a template bank, one can cover all masses in the 1-3 M Q range with some 
predetermined maximum loss in SNR [82, 83]. The smaller the maximum loss 
in SNR, the larger the number of templates needed in the bank. Typically, 
searches will implement a template bank with a maximum SNR loss of 3%, 
which leads to template banks containing of the order of a few hundred tem- 
plates (the exact number depends on the noise spectrum because both z(t) 
and a do, and therefore the number of templates can change from epoch to 
epoch). 
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When the noise is stationary and Gaussian, then matched filtering alone 
gives the best probability of detecting a signal (given a fixed false alarm rate) . 
However, as mentioned earlier, gravitational wave interferometer noise generi- 
cally contains noise bursts, or glitches, which provide a substantial noise back- 
ground for the detection of binary inspirals. It is possible for strong glitches 
to cause substantial portions of the template bank to simultaneously yield 
high SNR values. It is therefore highly desirable to have some other way of 
distinguishing the majority of glitches from true signals. 

The method which has become standard for this is to use a chi-squared 
(x 2 ) veto [84]. When a template exceeds the trigger threshold in SNR, it is 
then divided into p different frequency bands such that each band should yield 
1/p of the total SNR of the data if the high SNR event were a signal matching 
the template. The sum of the squares of the differences between the expected 
SNR and the actual SNR from each of the p bands, that is the \ 2 statistic, is 
then calculated. The advantage of using the x 2 veto is that glitches tend to 
produce large (low probability) x 2 values, and are therefore distinguishable 
from real signals. Thus, only those template matches with low enough \ 2 
values are considered triggers. 

If the data were a matching signal in Gaussian noise, the x 2 statistic would 
be x 2 distributed with 2p — 2 degrees of freedom [84]. However, it is much more 
likely that the template that produces the highest SNR will not be an exact 
match for the signal. In this case, denoting the fractional loss in SNR due to 
mismatch by /z, the statistic is distributed as a non-central chi-squared, with 
non-centrality parameter A < 2g 2 fi. This simply means that the x 2 threshold, 
X*, depends quadratically on the measured SNR, g 7 as well as linearly on \i. 

In practice, the number of bins, p, and the parameters which relate the % 2 
threshold to the SNR , as well as the SNR threshold g* which an event must 
exceed to be considered a trigger are determined empirically from a subset of 
the data, the playground data. A typical playground data set would be ^10% 
of the total data set, and would be chosen to be representative of the data set 
as a whole. Playground data is not used in the actual detection or upper limit 
analysis, since deriving search parameters from data which will be used in a 
statistical analysis can result in statistical bias. Values for these parameters 
for the LIGO SI and S2 BNS analyses are given in Table 1. 



Data Set 


8* 


P 


LI 


X* 


H1/H2 x* 


SI 

S2 


6.5 
6.0 


8 
15 


5(p + 
5(p + 


0.03£> 2 ) 
0.0l£> 2 ) 


5(p + 0.03f? 2 ) 
12.5 (p + 0.01 g 2 ) 



Table 1. Search algorithm parameters for SI and S2 BNS searches. These param- 
eters were determined using playground data extracted from the SI and S2 data 
sets respectively. Note that the \ 2 threshold, \* , is different for the Louisiana and 
Washington instruments in the S2 run. 
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Finally, let us say a few words about clustering. As discussed earlier, when 
a glitch occurs, many templates may give a high SNR. This would also be 
true for a strong enough signal. It would be a misinterpretation to suppose 
that there might be multiple independent and simultaneous signals - rather, 
it is preferable to treat the simultaneous events as a cluster and then try to 
determine the statistical significance of that cluster as a whole. The simplest 
strategy, and the one used thus far, is to take the highest SNR in the clus- 
ter and perform the \ 2 using the corresponding template. Another possibility 
might have been to take the template with the lowest x 2 value as representa- 
tive, or some function of g and \ 2 ■ In fact, there is reason to believe that the 
last option may be best [23]. 

3.3 Coincidence and Auxiliary Channel Veto 

Although the \ 2 veto reduces the rate of triggers from glitches, some glitch 
triggers survive. However, there are further tests that can be used to eliminate 
them. Most importantly, if more than one interferometer is involved in the 
search, one can require consistency between their triggers. LIGO is especially 
well designed in this regard. Because LIGO's three detectors are almost co- 
aligned, they should all be sensitive to the same signals (although the 2km 
H2 will only be half as sensitive to them as the 4km instruments). Thus, any 
signal that appears in one should appear in all. On the other hand, there is no 
reason for glitches in one instrument to be correlated with glitches in another, 
especially between a detector located in Washington state and the Louisiana 
instrument. Thus, one way to distinguish between triggers generated by actual 
gravitational waves and those generated by glitches is to demand coincidence 
between triggers in different instruments [85, 86]. 

The most fundamental coincidence is coincidence in time. The timing pre- 
cision of the matched filtering for properly conditioned intcrferometric data 
is ~ 1 ms. The larger effect is the time it could take the gravitational wave 
to traverse the distance between detectors (i.e. the light travel time between 
them). The maximum time delay for this is also measured in ms (e.g. 10ms 
between the Washington site and the Louisiana site for LIGO). Thus, trig- 
gers at one site which are not accompanied by triggers at another site within 
the light travel time plus 1 ms are likely not gravitational waves and can be 
discarded. For triggers from instruments which are sufficiently well aligned, 
there are several other quantities for which one could required coincidence. 
Of these, the only one which had been applied to date as a trigger veto is 
the template which generated the trigger - for the LIGO S2 search the same 
template was required to have generated all coincident triggers (or represent 
all coincident clusters of triggers) or they were discarded. 

The final hurdle that a trigger may have to overcome to remain viable 
is that it not be associated with a known instrumental disturbance. Auxil- 
iary channels which monitor the instruments and their environment contain 
information about many potential disturbances. Those channels most likely 
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to correspond to spurious disturbances which would be manifest in the grav- 
itational wave channel have been studied intensively. To date, studies have 
revealed that channels which allows for safe and useful auxiliary channel ve- 
toes are not forthcoming for most instruments. However, for LIGO's second 
science run it was discovered that a channel which measures length fluctua- 
tions in a certain optical cavity of the LI interferometer had glitches which 
were highly correlated with glitches in that instrument's gravitational wave 
channel [87]. Thus, triggers which occurred within a time window 4 seconds 
earlier to 8 seconds later than a glitch in this auxiliary channel in Louisiana 
were also discarded for that analysis. 

Finally, triggers which survive all of these cuts must be examined indi- 
vidually to determine if they are genuine candidates for gravitational wave 
signals. The flow chart for the procedure we have just describe is shown in 
Fig. 5. A more comprehensive and detailed discussion of such an algorithm 
can be found in [88]. Details of the pipelines used by LIGO for the SI and 
S2 analyses can be found in [89] and [90] respectively. In the next section, we 
discuss how the surviving candidates are analyzed and how upper limits are 
set from them. 



4 Statistics and Results 
4.1 Background 

Once the method outlined above is completed, one is left with zero or more 
coincident triggers that need to be interpreted. The pipeline, and in particular 
the search parameters g* and x* , are chosen so that the probability of missing 
a real signal are minimized. This means, in practice, that the probability of 
having coincident triggers from glitches is not minimized. 

These coincident triggers from glitches form a background against which 
one does a statistical analysis to determine the likelihood that there are signals 
(foreground events). In order to accomplish this task, it is useful to understand 
the background rate. Ideally, one would have a model for the background, but 
there is no such model for GW detector glitches. However, there is a fun- 
damental difference between background and foreground triggers that can 
be exploited - coincident background triggers are caused by random glitches 
which happened to coincide between detectors (assuming that the glitches 
in detectors are not correlated), where as foreground triggers are coincident 
because they are caused by the same gravitational wave signal. Thus, if one 
introduces a large enough artificial time delay between detector data streams 
(e.g. greater than 11ms between the triggers from Louisiana and those from 
Washington), the background rate should remain unchanged while the fore- 
ground rate must vanish. 

A large number of time shifts (e.g. 40 for S2) are used to provide a great 
deal of background data and thus increase confidence in the inferred back- 
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Fig. 5. A diagram showing a typical work-flow to look for signals from neutron 
star binaries. This is a two instrument work-flow for LIGO involving one detector 
from Washington and the Louisiana detector. It is similar to the one used for the S2 
analysis. There would a slightly more complicated work-flow diagram when all three 
instruments were used. Note that the Washington data is only analyzed at times and 
for binary templates that correspond to a trigger in Louisiana, thus minimizing both 
false alarms and computing time. The + and - signs by the two bottom-most arrows 
indicated that the addition and subtraction of coincident triggers respectively - i.e. 
passing the SNR and \ 2 thresholds adds coincident triggers, while occurring during 
an auxiliary channel glitch removes them. 



ground rate. Both the mean and variance of the background rate are esti- 
mated. The foreground rate is then compared. If the foreground rate is signif- 
icantly higher than the mean background rate, it is possible that some of the 
coincident triggers are from gravitational waves. A careful look at each coinci- 
dent trigger then ensues. Criteria such as consistency of SNR values, problems 
with auxiliary channels, and behaviour of the gravitational wave channel are 
invoked to try to eliminate surviving background events. In particular, it is 
known that poor instrument behavior often leads to many consecutive glitches, 
and coincident triggers are much more likely during these glitchy times. 
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4.2 Upper Limits and the Loudest Event 

If no gravitational wave candidate is identified amongst the coincident trig- 
gers, then the best one can do is set upper limits on the rate at which such 
signals occur. Clearly, it is meaningless to quote a rate unless one specifies the 
minimum SNR one is considering. Since triggers can have different SNR's in 
each instrument, for coincident searches it is convenient to consider a measure 
of the combined SNR. 

If the noise in all instruments were Gaussian, then it would be appropri- 
ate to combine SNR's in quadrature, e.g. g c = g 2 H + q\, where Qh and Ql 
are the SNR's at Washington and Louisiana respectively. This would still be 
appropriate for non-Gaussian noise if the rate and SNRs of glitches were the 
same for instruments at both sites. However, this is not guaranteed to be the 
case, and if it is not, the noisier instrument can dominate the the upper limit 
statistics. To avoid this problem, the combined SNR needs to be modified 
slightly to 

Q 2 c = Q 2 H +a Q l, (17) 

where a needs to be determined empirically. 

Now, let us denote the mean rate at which signals arrive at a set of de- 
tectors with SNR qc > Q* c by 1Z. If we model the arrival of such signals as a 
Poisson process, then the probability of detecting such a signal within time T 
is given by 

P{QOQ*ci H) = l-e nT <e*c\ (18) 

where e(g* c ) is the detection efficiency, i.e. the ratio of detected signals to 
incident signals at the threshold q* c . This is not the probability of observing 
a trigger with qq > Qc however, because it does not account for background 
triggers. If the probability of having no background event with qc > Q* c 1S Pb, 
then the probability of observing at least one trigger with qc > Q* c is 

P( 9C > g c ;ll,b) = l- P b e nT ^h). (19) 

The loudest event (i.e. the event with the maximum combined SNR, g max ) 
sets the scale for the upper limit on the rate. More precisely, to find the 90% 
frequentist upper limit, one needs to determine the value of 1Z such that there 
is a 90% chance that no combined SNR would exceed p max over the course of 
the run. In other words, one needs to solve 0.9 = P(qc > ftnax ; H, b) for 1Z. 
Doing so, we find 

_ 2.303 + lnA . . 

'M30% - \ * (20) 

There are two quantities in (20) which need to be ascertained. The first is the 
efficiency of the detection algorithm to signals with combined SNR g max . This 
is evaluated through Monte Carlo simulations where simulated signals from 
the target population of neutron star binaries (using the population models 
discussed in Sect. 2.3) are injected into the data. 
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The second is the background probability, P b . The most straightforward 
approach might be to estimate Pb using the background events resulting from 
the timeshifts as described in Sect. 4.1. However, these background rates are 
known to be subject to significant variation depending on the details of search 
- e.g. reasonable changes to the event clustering criteria lead to different 
background rate estimates [23]. Since for detections one would follow up with 
detailed investigation of coincident triggers anyway, the background estimates 
are used in a more qualitative manner, and this variation is not an issue for 
detection. For an upper limit, however, one needs a quantitative result for the 
background rates, and if there is uncertainty in the value obtained, it also 
must be quantified. Failing to do so could lead to undercoverage, i.e. a 90% 
upper limit which is below the actual rate that can be inferred from the data. 
Undercoverage is considered a "cardinal sin" in frcqucntist analyses. 

Therefore, rather than try to get a quantitative estimate, the standard 
practice is to simply use Pb = 1 in (20). Note that this maximizes 1Zqo% with 
respect to Pb- It therefore gives a conservative upper limit - the actual 90% 
confidence rate is certainly lower. While undesirable, this is considered a "ve- 
nial sin" in frequentist statistics, and therefore far preferable to undercoverage. 
Furthermore, since Pb is the probability that no background coincidences will 
occur with SNR above g max , and since g max is the SNR of the loudest back- 
ground coincidence that actually did occur, it is statistically unlikely that the 
actual TZqo% would be very much below the one obtained by this method. 
This statistic is known as the loudest event statistic, and is discussed at some 
length in [91]. 

Finally, BNS rate limits are typically quoted in units of "per Milky Way 
equivalent galaxy (MWEG) per year". This results from expressing the effi- 
ciency in terms of effective number of Milky Way equivalent galaxies to which 
the search was sensitive, Nq- The conversion between e(g max ) and Nq is 



where L pop is the effective blue-light luminosity of the target population and 
Lq = 9x 10 9 Lq is the blue-light luminosity of the Milky Way galaxy. In terms 
of No, the 90% frequentist rate upper limit is written 




(21) 



K 90% = 2.303 x ( ±£) x yr" 1 MWEG 




-l 



(22) 



where T has units of years and Pb = 1 has been used. This, then, is the most 
common form of the upper limit that is quoted in the literature, and in this 
article. 



Searches for Gravitational Waves from Binary Neutron Stars: A Review 23 
4.3 Results 

To date, there have been five published searches for BNS coalescence by large 
scale (> 100 m) interferometric detectors 6 . The first was performed by TAMA 
with a single instrument [26] on the 1999 DT2 (Data Taking 2) data. Subse- 
quent multi-detector searches were performed by LIGO based on SI [25] and 
S2 [23] data and jointly by LIGO and TAMA [20]. Most recently, TAMA has 
analyzed cumulatively data taken over many Data Taking runs between 2000 
and 2004. No detections have been claimed for any analysis thus far. The re- 
sults therefore have all been observational upper limits on coalescence rates. 
They are summarized in Table 2. 



Data Set 


T 


N G 


^•90% 




(hrs) 




(MWEG-iyr- 1 ) 


TAMA DT2 


6 




5170 


LIGO SI 


236 


0.6012:13 


170 


LIGO S2 


339 


l-34±g;g? 


47 


LIGO S2/TAMA DT8 


584 


0.761°;°° 


49 


TAMA 2000-2004 


2075 




20 



Table 2. Results from five searches for neutron star binaries involving data from 
large scale (> 100 m) interferometric gravitational wave detectors. T is the obser- 
vation time used for setting the upper limit, i.e. the total observation time minus 
playground data set length minus time periods lost due to vetoes. In the first and 
last entry, TAMA did not convert their efficiency to Ng- The results quoted in these 
rows, therefore, are not in units of MWEG -1 yr _1 , but rather in "per population 
observed yr _1 ". Note that where Ng is used, the quoted rates are those from the 
lower bound for Ng, thus giving the most conservative upper limit. 



All of these rates are substantially above the theoretically predicted rates 
for these searches [92]. Nonetheless, these are the best direct observational 
limits on neutron star binaries to date. Furthermore, even at these sensitivities, 
given how little is known about these gravitational wave sources, there is a real 
chance, however marginal, of a serendipitous discovery. Finally, these studies 
were used as a testing ground to develop the analysis tools and methodologies 
that will maximize the possibility of discoveries in future analyses. 



6 One other was done using TAMA300 and LISM, a 20 m prototype 
interferometer [24] - we will not consider this one here since we feel that the 
results are essentially subsumed by the TAMA 2000-2004 result 
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5 Future Prospects 

5.1 Interferometers Now and Future 

As mentioned in the introduction, to date (Spring 2007), LIGO has completed 
four science runs, SI from 23 Aug-9 Sep 2002, S2 from 14 Feb-14 Apr 2003, 
S3 from 31 Oct 2003-9 Jan 2004, and S4 22 Feb- 23 Mar 2005. These short 
science runs were interspersed with periods of commissioning and have yielded 
dramatic improvements in sensitivity (see Fig. 3). LIGO is now operating with 
its design sensitivity. The fifth science run, S5, began in Nov 2005, and has the 
goal of collecting a full- year of coincident data at the LIGO design sensitivity. 
This goal is expected to be achieved by late summer 2007 [62]. The scientific 
reach of S5 - in terms of the product of the volume of the Universe surveyed 
times the duration of the data sample - will be more than two orders of 
magnitude greater than the previous searches. 

A period of commissioning following the S5 run will hopefully improve 
the sensitivity of LIGO by a factor of ~ 2 [62, 61]. An anticipated S6 run 
with the goal of collecting a year of data with this "enhanced" LIGO interfer- 
ometer is projected to commence in 2009. If S6 sensitivity is indeed doubled, 
enhanced LIGO will survey a volume eight times as great as the current LIGO 
sensitivity. 

Following S6, in 2011, the LIGO interferometers are slated for decommis- 
sioning in order to install advanced interferometers [62, 93]. These advanced 
LIGO interferometers are expected to operate with 10 times the current sen- 
sitivity (a factor of 1000 increase in the volume of the Universe surveyed) by 
- 2014. 

In addition to LIGO, the Virgo and GEO600 observatories are also par- 
ticipating in the S5 science run and are expected to undergo upgrades that 
will improve their sensitivity along with LIGO. In Japan, TAMA will hope- 
fully be replaced by the Large Scale Cryogenic Gravitational-wave Telescope 
(LCGT), a subterranean cryogenic observatory with similar performance to 
advanced LIGO [62]. Having several detectors operating at the time of a grav- 
itational wave event will allow better determination of an observed system's 
parameters. 

5.2 Future Reach and Expected Rates 

It is possible, even before doing analyses, to get estimates of how instruments 
will perform in future analyses given their noise curves. For instance, since S2 
a sequence of improvements to LIGO's sensitivity have been made (see Fig. 3) 
that have resulted in an order of magnitude increase in the range to which 
it is sensitive to neutron star binary inspiral. Galaxy catalogs can be used 
to enumerate the nearby galaxies in which BNS inspirals could be detected. 
The relative contribution to the overall rate of BNS inspirals for each galaxy 
is determined by its relative blue light luminosity compared to that of the 
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Milky Way. Monte Carlo methods are used to determine the fraction of BNS 
signals that would be detectable from each galaxy during a particular science 
run. As described above, this procedure yields a figure for sensitivity during 
any science run: the number of Milky Way Equivalent Galaxies (MWEG) that 
are visible. 

In the case of the second LIGO science run, S2, the search was sensitive to 
1.34 MWEG. The S2 run produced a total of 339 hours (around 0.04 years) 
of analyzed data. The inverse of the product of the number of galaxies to 
which a search was sensitive times the livetime is a measure of the scientific 
reach of the search. For S2 this was ~2x 10 7 Myr" 1 MWEG" 1 . For initial 
LIGO sensitivity (which has been achieved during the current S5 science run) 
Nutzman et al. [94] estimate the effective number of MWEG surveyed to be 
—600 MWEG, though this number depends on the actual sensitivity achieved 
during the search. Assuming that S5 produces one year of analyzed data, 
the scientific reach of this search is —2000 Myr" 1 MWEG" 1 . These numbers 
should be compared to estimated BNS merger rates in the Milky Way. For 
example, the "reference model" (model 6) of Kalogera et al. [92] quotes a 
Galactic rate of between -20 and -300 Myr -1 MWEG" 1 with a most likely 
rate of -80 Myr -1 MWEG" 1 . 

Extrapolating to the future, the enhancements that are expected to double 
the range of LIGO before the sixth science run could increase the scientific 
reach of LIGO to -200 Myr" 1 MWEG" 1 , which would now begin to probe 
the expected range of BNS coalescence rates. Since the Advanced LIGO is ex- 
pected to improve the sensitivity by a factor of 10 compared to the current sen- 
sitivity, given a few years of operation a scientific reach of — 1 Myr" 1 MWEG" 1 
should be achieved. It is therefore likely that detection of BNS inspirals will 
become routine during the operation of Advanced LIGO. 

5.3 BNS Astrophysics with GWs 

As discussed above, we expect that future GW observations will begin to probe 
the interesting range of BNS inspiral rates over the next several years; when 
advanced detectors like Advanced LIGO begin running at the anticipated 
sensitivity, we expect BNS inspirals to be routinely detected, which will give a 
direct measurement of the true BNS merger rate as well as observed properties 
(such as the distribution of masses of the companions) of the population. Such 
constraints on the population of BNS can then be compared to the predictions 
from population synthesis models, which then can give insight into various 
aspects of the evolution of binary stars [71]. 

Design differences between current and future interferometers might also 
open new avenues of investigation. One intriguing possibility involves mak- 
ing measurements of neutron star equations of state using GW observations 
of BNS mergers. It is possible that advanced interferometers will be able to 
be tuned to optimize the sensitivity for a particular frequency band (while 
sacrificing sensitivity outside that band). This is already being discussed as 
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a feature for Advanced LIGO. This raises the possibility of tuning one of the 
LIGO interferometers to be most sensitive to gravitational waves at high fre- 
quencies, around 1 kHz, where effects due to the size of the neutron stars is 
expected to be imprinted in the gravitational waveform. Advanced interfer- 
ometers could then make a direct measurement of the ratio of neutron star 
mass to radius and thereby constrain the possible neutron star equations of 
state [95]. 

Finally, recent evidence suggests that short hard gamma-ray bursts (GRBs) 
could be associated with BNS mergers or neutron-star/black-hole binaries. 
Once GW interferometers are sufficiently sensitive, a search for inspiral wave- 
forms in conjunction with observed GRBs could confirm or refute the role of 
binaries as GRB progenitors if a nearby short GRB were identified. If the bi- 
nary progenitor model is accepted then gravitational wave observations of the 
associated inspiral could give independent estimates or limits on the distance 
to the GRB. While most GRBs will be too distant for us to hope to detect 
an associated inspiral, a few may occur within the observed range. 

6 Concluding Remarks 

We have reviewed here the recently published results of searches for grav- 
itational waves from neutron star binaries. No detections have been made 
in analyses published thus far, but this is hardly surprising given the gap 
between current sensitivities and those required to reach astrophysically pre- 
dicted rates. Nonetheless, steady progress is being made in refining instru- 
mentation and analyses in preparation for that time in the not-too-distant 
future when this gap has been closed. 

ft should be obvious to the reader that this has hardly been a comprehen- 
sive description of any aspect of this search. Indeed, to provide such a descrip- 
tion would require at least an entire volume in itself. Notably, our description 
of the mathematical theory of gravitational waves, the theoretical and com- 
putational underpinnings of population synthesis modeling and our discussion 
of interferometric design were all sketchy at best. Nor have we delved at all 
into the important topics of error analysis and pipeline validation. 

Nonetheless, we have attempted to provide a birds-eye view of the relevant 
aspects of searches for gravitational waves from coalescing neutron star bina- 
ries with enough references to the literature to provide a starting point for 
readers interested in any particular aspect. And there is good reason to believe 
that, in time, gravitational wave searches will become increasingly interest- 
ing and relevant to a ever broadening group of astronomers and astrophysi- 
cists outside the gravitational wave community. Already, viable (if somewhat 
marginal) theoretical models have been constrained [10], and greater volumes 
of more sensitive data are in hand. It seems that it is only a matter of time 
before GW detectors begins seeing the Universe through gravitational waves. 
We look forward to the opportunity to review the findings when they do. 
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